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ABSTRACT 

Energy versus magnetic- field diagram of the spin-1 Haldane system 
with an impurity bond is studied in terms of spin- 1/2 degree of free- 
dom at the sites neighboring the impurity bond by means of analytical 
method. We examine the equivalence between the realistic Hamiltonian 
and the phenomenological Hamiltonian which is composed two spin- 1/2 
spins representing the spin- 1/2 degree of freedom. It is proved that when 
the strength of the impurity bond is sufficiently weak, the two Hamiltonians 
are equivalent to each other, as far as the energies of the low-lying states 
are concerned. We determine the correspondence between the interaction 
constants in the phenomenological Hamiltonian and those in the realistic 
Hamiltonian. 
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The impurity effects on the quantum antiferromagnetic chain have been 
subject of a large number of theoretical and experimental studies. One of 
the recent exciting topics of this subject is the edge and/or impurity ef- 
fects on the spin-1 antiferromagnetic Heisenberg chain which is the simplest 
system realizing Haldane's conjecture[l] of the difference between integer- 
spin and half- integer-spin antiferromagnetic Heisenberg chains. Kennedy [2] 
found that the open chain has a fourfold degenerate ground state composed 
of a singlet and a triplet which we call the Kennedy triplet, in contrast to a 
unique singlet ground state of the periodic chain [3]. The fourfold degener- 
acy of the ground state, which was originally found in the so-called AKLT 
model [3] with open boundary conditions, is considered to reflect the hidden 
Z2XZ2 symmetry in the open chain [4]. The hidden symmetry is attributed 
to the spin- 1/2 degree of freedoms at edges of the chain. Recently the 
present authors and Harada[5] investigated theoretically the the impurity- 
bond effect on the ground state and low-lying excited states of the chain 
in terms of domain-wall excitations to interpolate between the open-chain 
and the periodic-chain cases. They showed that the impurity bond brings 
about a massive triplet mode in the Haldane gap and that the triplet state 
comprises three of the four ground states of the open chain. Miyashita and 
Yamamoto[6] perform the Monte Carlo analysis of the open chain to show 
that the magnetic moments localized around the edges for the Kennedy 
triplet decay exponentially with the decay constant which is about 6 in 
lattice spacing. On the experimental side, on the other hand, Kikuchi et 
a/. [7] investigated the Zn impurity effects on the magnetization process of 
the spin-1 Haldane system NiC2042DMIz. Hagiwara et a/. [8] performed 
the ESR experiment on Ni(C2H8N2)2N02(C104), abbreviated NENP, con- 
taining a small amount of spin- 1/2 Cu^+ impurities and gave experimental 
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evidence for the existence of the spin- 1/2 degrees of freedom at the sites 
neighboring the impurity spin. They analyzed successfuhy their experimen- 
tal results by using the phenomenological Hamiltonian composed of three 
spin- 1/2 spins with the anisotropic exchange interactions. The present au- 
thors very recently proved the equivalence between the phenomenological 
Hamiltonian and a more realistic Hamiltonian and gave a clear interpreta- 
tion of the anisotropy of the exchange interaction [9]. 

The purpose of this article is to investigate the energy versus magnetic- 
field diagrams of the spin-1 Haldane system with an impurity bond in 
terms of the spin- 1/2 degree of freedom by means of analytical method. 
We express the Hamiltonian Ti of the chain with an impurity bond under 
the magnetic field as sum of TYq, 'H' and T^mf as 

Ti — Tio + + Hmf 5 (la) 

no=J2 + , (lb) 

i=l i=l 
H' = HihN,i , (Ic) 

h^^e = J[SfSf, + Spl + XSISI.) (J, A > 0) , (Id) 

^mf = -p(^a;5't^otal + ^^'S'total) » (1^) 

where Si (^ = 1, 2, • • •, N) is the spin-1 operator; ^^o^a-i = ^(_=iSf and 
'^fotai = Yl!i=i^^'^ 9 is ^-factor; A is the parameter which represents 
the anisotropy of the exchange; k, represents the strength of the impurity 
bond, with At = and k= \ corresponding, respectively, to the open chain 
and the periodic chain. Thus, Tio, H' and Hmi represent, respectively, the 
Hamiltonian for the coupling through the host bonds, that for the coupling 
through the impurity bond, and that for the Zeeman energy. 
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In ref. 5, we calculated the energies of the low- lying states of Ho-\-H' 
by the variational method, assuming that |/^| <C 1, that is, treating ?{' as 
a small perturbation. Using this method we calculate the energies of the 
low-lying states of H. In the Haldane region, to which we confine ourselves 
hereafter, the ground state of the unperturbed Hamiltonian Tio is fourfold 
degenerate [2, 3]. Since the fourfold degeneracy reflects the existence of the 
spin- 1/2 degree of freedom at the neighboring sites of the impurity bond, we 
may expect, as in the previous paper [9], that the energies of the low-lying 
states of the present system can be calculated from the phenomenological 
Hamiltonian given as 



where si and S2 are the spin operators which represent the spin- 1/2 degrees 
of freedom at the sites neighboring the impurity; J is the effective exchange 
constant; A is the parameter which represents the anisotropy of the effective 
exchange; g^ and gz are the effective ^f-factors. By comparing the energy 



-^phe equivalent to 7i, as far as the energies of the low-lying states are 
concerned. 

The variational method discussed in ref. 5, according to which the Hal- 
dane region is given by 4 > d > 2A — 4, leads to the following fourfold 
ground-state wave functions $ and ^^^^ (r = -|-, 0, — ) expressed in the 
matrix-product form[10,ll]. The function $ describes the state with no 
domain wall and is given by 



expressions for T-^P^e with those for H, we will show that, when \k\ « 1 



$ = Trace [0i02 • • • 0Ar-i0iv] 

A sin(6' 
(pi = cos (6') Ci az H 




(3) 



(4) 
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where ai, Q, and f3£ are the spin states at the ^-th site, which correspond, 
respectively, to 51 = 1 , 0, and — 1 , and a± [= (ax^idy ) / \/2 ] , ax, CTy, and cr^ 
are the Pauh matrices. The parameter 9 is determined from the equation, 

cos(2e) = • (5) 

(t) 

The function describes the states with a domain wall and is expressed 
in terms of and the wall operator w as 

= Trace [0102- --^AT- 107V w;] , (6) 

where w = —a- for r = +, w = cr^ for r = 0, and w = (j+ for r = — . It is noted 
that, when A = l and = 0, $ represents the singlet state, and ^n^, 
and represent, respectively, the triplet (Kennedy triplet) states with 
S^Q^^i = l, 0, and —1. It should be noted that a Neel-type configuration of 

(t) 

magnetic moment near the edges is accompanied by the state • ^^^^ 

energy expectation value Es{k) ( E^'^'' ) of Ho+H' with respect to the state 
$ ( ^^jp ) is given as 

EM = -iN-l + .)J{ (1^) (l + J) } + NM^^ . (7) 



Ei"^ = E,{k) + Ao{k) 

(4- 

(4 -A) 



= E,{k) + «J ^-^^ >-} , (8) 

„/x ^ , (i - d){8 + 2d - \d) -, 
= E,(k) + kJ{^- '-} . (9) 

We now calculate the energy versus magnetic field diagram under the 
same assumption |k| <C 1 as previous. By using expression of the matrix 



elements given by 



(^l^totl^) = (^i?Vtotl^S?^) = 0, (10a) 
(^isT^l^totl^Sv^^) = -i^^N^lSLtl^^N^) = 1 , (10b) 

= + = (12 + .-4A) - 



(^!v^^l^fotail^^"^) = 0, (lOe) 

we obtain the secular equation for the energy difference (the energy 
measured from Es{k) ) as 



AE 



[Ao{k) - AE] [{A^{k) - AeY - [gH.f] 

-2{GgH,f{Ai{ti) - AE]] = . (11) 



The energy difference AE versus magnetic field Hx or Hz diagram within 
this approximation can be constructed by solving this equation [12]. 

Let us examine the relation between Ti^^^ and 7i. The four energy 
eigenvalues and E^'^^ (r = 0, ±) of H^^^ for 11^ = and = are 
calculated to be 



=-^(2 + A), (12) 



=E, + Ao = E, + J, (13) 

^i^) =E, + Ai = Es^ ^(1 + A) , (14) 

where the subscripts s, t and the superscript r have the same meaning as 
those of eqs. (7) to (9) but for the two-spin-1/2 system. The equation for 
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the energy difference /S.E (the energy measured from ) is easily obtained 

as 

A^[[Ao - A^;] [{Ai - A^}' - [gAf] 

-2(^,i?,)2{Ai-A^}] = 0. (15) 

Comparing eqs. (12)-(14) and (15) with eqs. (7)-(9) and (11), we can de- 
termine the correspondence between the interaction constants in HP^^ and 
those in 7i. The results are 



2{A-d){A-2\ + d) 



"- 2(4-2A + c/) - ^^'^ 

9z^ 9 (18) 
„ V2(4 - d)(4 - 2A + 4 

^ (12 + d-4A) ^ ^ 

We have thus shown that, TC^^^ is equivalent to 7i when |/«| <C 1, as far as 
the energies of the low- lying excited states are concerned. Equations (17)- 
(19) give a clear explanation for the origin of the anisotropy in Ti^^'^; as 
seen from eqs. (17)- (19), in the case of A = A' = 1 the uniaxial anisotropy d in 
H produces the anisotropy in H^^^. This result confirms again the concept 
of the spin-1/2 degrees of freedom at edges of the open spin-1 chain. The 
use of W^^^ for the semi-quantitative analysis of the experimental results 
is also justified. 

Finally we discuss the solutions of eq. (11). The equation has a trivial 
solution 



AE = AEs = 0. 
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(20) 



The three solutions other than this trivial solution , which we denote as 
AE^'^\k) (r = 0, zb), are obtained for the case of Hx = and Hz=H >0 
as 

AeI''\k) = Ao{k) , (21) 
AE^^^ (/^) = Ai + rg^H, , (r = ±) . (22) 

Those for the case of Hx > and Hz = are obtained as 

Asf (/t) = Ai(k) , (23) 



r^i^^^^^^^^y + iGgH^y (r = ±). (24) 

In Fig. 1 we show the AE versus magnetic-field diagram calculated from 
either eqs. (21)-(22), eqs. (23)-(24), or eq. (11) for the case of d>0 and k>0. 
It should be noted that we can calculate the magnetic-field or temperature 
dependence of the Neel-type configuration of magnetic moment around 
impurity bond at low temperatures and that the results obtained in the 
present paper are in qualitatively good agreement with those obtained by 
the numerical diagonalization. Details of the calculations will be published 
in the near future. 
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Figure Caption 

. 1. Schematic energy difference /S.E (the energy measured from E^) 
versus magnetic-field diagram for the case of d > and > obtained 
by the present analysis; (a) for = and = H > 0, (b) for 
H^ = H>and H^ = 0, (c) for H^ = H^ = ^>0. 
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